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Abstract 

We examine and present new combinatorics for the Schur polynomials from the view¬ 
point of quantum integrability. We introduce and analyze an integrable six-vertex model 
which can be viewed as a certain degeneration model from a t-deformed boson model. 
By a detailed analysis of the wavefunction from the quantum inverse scattering method, 
we present a novel combinatorial formula which expresses the Schur polynomials by us¬ 
ing an additional parameter, which is in the same sense but different from the Tokuyama 
formula. We also give an algebraic analytic proof for the Cauchy identity and make appli¬ 
cations of the domain wall boundary partition functions to the enumeration of alternating 
sign matrices. 


1 Introduction 

Schur polynomials is an ubiquitous object in mathematics and mathematical physics, rang¬ 
ing from representation theory, combinatorics, enumerative geometry to knot theory. Being 
one of the most fundamental symmetric polynomials, extensive studies and numerous com¬ 
binatorial identities have been found for the Schur polynomials. Many mathematical objects 
are introduced for which the Schur polynomials serves as a building block since one can use 
various properties found for the Schur polynomials. One typical example is a stochastic pro¬ 
cess called the Schur process [Tj whose probability measure is given by the Schur polynomials, 
from its property many useful combinatorial formulae can be employed to study the dynamics 
of the process. 

The Schur polynomials itself has many expressions. One of the most famous ones is the 
Jacobi-Trudi identity which expresses Schur polynomials in terms of elementary symmet¬ 
ric polynomials. Besides these traditional identities, a very interesting formula called the 


* E-mail: motegi@gokutan.c.u-tokyo. ac.jp 
'E-mail: sakai@gokutan.c.u-tokyo. ac.jp 


1 



Tokuyama formula Emm was found which expresses Schur polynomials in a combinatorial 
form. The distinctive feature of the formula is that the expression is given in terms not only 
of the symmetric variables but also with an additional parameter which does not appear in 
the original definition of the Schur polynomials. Specializing the additional parameter, the 
Tokuyama formula reduces to the determinant formula, Weyl character formula and so on. 
In this sense, the Tokuyama formula can be viewed as a deformation of Weyl character. 

Recently, the Tokuyama formula has found its interpretation in terms of integrable vertex 
models [5]. The idea to understand this Tokuyama formula from the point of view of quan¬ 
tum integrability is to extend the L-operator for the free-fermion model (the corresponding 
R-matrix is the trigonometric Felderhof model) to an L-operator including an additional 
parameter besides the spectral parameter. The newly-introduced parameter plays the role 
of refining the combinatorial expression for the Schur polynomials. This idea of relating 
Tokuyama formula which is a deformation of Weyl character formula to integrable vertex 
models have been recently generalized to the factorial Schur polynomials [6] and other types 
of symmetric polynomials mmmm by introducing inhomogeneous parameters or changing 
boundary conditions. 

In this paper, we present another type of combinatorial formula for the Schur polynomials 
with an additional parameter by investigating a partition function of an integrable six-vertex 
model. We find the following combinatorial formula for the Schur polynomials 
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where /3 is an arbitrary parameter, x^ = ( x ..., x^), k = 1,... ,N are strict partitions 
satisfying the interlacing relations x^ >- x^ N ~ 1 ^ >-■■■>- x^ = 4>, x^ N ^ is fixed by the Young 
diagram A = (Ai,..., \n) as A j = x^ — N + j — 1, and #(y\x) denotes the number of parts 
in y which are not in x. 

We obtain this combinatorial formula from an integrable six-vertex model. The six-vertex 
model we consider can be regarded as a certain degeneration of a f-deformed non-Hermitian 
boson model introduced in m- The corresponding wavefunction seems to be analyzed by 
the coordinate Bethe ansatz [12j which at the degeneration point is given by the Schur 
polynomials times a factor including an additional parameter. To obtain the combinatorial 
formula (II.ID . we investigate the corresponding model from the point of view of the quantum 
inverse scattering method, i.e., starting from the L-operator which is the most fundamental 
object in quantum integrable models. We view the wavefunction as a partition function of 
a six-vertex model and evaluate it in two ways. First, we evaluate the partition function 
directly to show that it is given by a Schur polynomials times an additional factor. Another 
way of evaluation is to take the viewpoint that the partition function consists of a layer of 
L-operators, and calculate the matrix elements of a single B-operator to show a summation 
formula for the partition function. Combining the two expressions give the combinatorial 
formula (11.11) . 

Besides the combinatorial formula, we examine the scalar products and present an alge¬ 
braic analytic proof for its determinant form, which combined with the wavefunction gives 
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us the celebrated Cauchy formula for the Schur polynomials. We also make an application 
of the domain wall boundary partition function to show a simple formula for a particular 
generating function of alternating sign matrices. 

This paper is organized as follows. We introduce the integrable six-vertex model as a 
particular degeneration of a non-Hermitian f-deformed boson model in the next section. We 
give an algebraic analytic proof for the scalar product in section 3. In sections 4 and 5, we 
evaluate the matrix elements and the wavefunction of the six-vertex model. We combine the 
results of sections 3, 4 and 5 to present the combinatorial formulae for the Schur polynomials 
in section 6. In section 7, we give an application of the domain wall boundary partition 
function to a generating function of the alternating sign matrices. 

2 Non-Hermitian ^-deformed boson model and reduction to 
the six-vertex model 

We introduce the f-deformed non-Hermitian boson model mm and its equivalent six-vertex 
model at the point t = — 1. The model is characterized by f-deformed boson algebra with 
generators N, B,B^ satisfying the following relations 

[B, = t N (l — f), [N,B] = -B, [N,B'] = Bl (2.1) 

These generators act on a Fock space T spanned by the orthonormal basis | n) (n > 0) as 

B\n) = |n — 1), B^\n) = (1 — t n+1 )|n + 1), N\n) = n|n). (2.2) 

The operators the dual orthonormal basis (n| (n > 0) satisfying (n|m) = 5 nm as 

(n|L^ = (1 — f n )(n — 11, (n\B = (n + l\, (n\N = n(n\. (2-3) 

We consider the following L-operator which is a slightly deformation of the one in m 

(2.4) 

acting on the tensor product W a ® J-j of the complex two-dimensional space W a and the Fock 
space at the jth site Tj. Let us denote the orthonormal basis of W a and its dual as (|0) a , |l) a } 
and {a(0|,a(l|}- The explicit forms of the matrix elements of (12.411 in the orthonormal basis 
are 


£aj(v) = 


1 - vB 

B,; v 


a(0\j{mj\C a j(v)\0) a \nj}j = (1 - /3vt n i)5 mjtnj , (2.5) 

a (1 \j \£ , ajiy) |0)a — &mj,nj— 1> (2-6) 

a(01j (rrij | C a j(v)\l) a \rij)j = v(l - t nj+1 )6 mjin:j+ i, (2.7) 

a (1 Ij (?7lj |71 a j (n) | l) a — v6• (2-8) 

In figure [TJ we depict the non-zero elements of the L-operator. 

The L-operator satisfies the intertwining relation (LLL-relation) 

{y^B a j (u)L a ^('u, u), (2.9) 
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Figure 1: The non-zero elements of the L-operator. 


which acts on W a <g> W b <g> ■ The L-matrix 
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is a solution to the Yang-Baxter equation: 

R ab (u,v)R ac (u,w)R bc (v,w ) = R bc (v,w)R ac (u,w)R ab (u,v). 


( 2 . 10 ) 


( 2 . 11 ) 


From the L-operator, we construct the monodromy matrix 

M 

Ta(v) = Y[C aj (v), (2.12) 

1=1 

which acts on W a <S> {R\ <8> • • • <8 'iFm)- Tracing out the auxiliary space, one defines the transfer 
matrix t(v) G End(J r ® M ): 


t(v) = Tv Wa T a {v). 


(2.13) 


Simplifications happen for the f-deformed boson model at t = — 1. Due to the matrix 
element flZZZD , when we construct A^-particle states from the vacuum, each site can only be 
occupied by at most one boson. This is equivalent to the reduction of the f-deformed boson 
model to the six-vertex model whose L-operator is given by 


£ a j(y) 


fl-/3v 0 0 0\ 

0 1 + f3v 2v 0 

0 1 v 0 
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(2.14) 


One can check that (12.141) indeed satisfies the intertwining relation (12.91) . We analyze the 
structure of this six-vertex model in this paper. We remark that at the other degenerated 
point f = 0 of the f-deformed boson model, the model is called the non-Hermitian phase model 
m , whose wavefunction m constructed from the L-operator is given by the Grothendieck 
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Figure 2: The diagrammatic representation of the elements of the monodromy matrix (13.11) 
with inhomogeneous parameters WM+i-j associated with site j. 


polynomials na nsi nzi ns hsi, which furthermore reduces to the Schur polynomials by 
taking f3 = 0 where the corresponding model is the Hermitian phase model mmmm- 
On the other hand, the wavefunction at [3 = 0 for generic t is given by the Hall-Littlewood 
polynomials mm- 


3 Scalar Products of state vectors of the six-vertex model 


In this section, we first construct a state vector of the integrable f-deformed boson model 
following the standard procedure of the quantum inverse scattering method (i.e. the algebraic 
Bethe ansatz) which is based on the Yang-Baxter algebra. The IV-particle state |({^}at)) 
is characterized by N unknown numbers Vj G C (1 < j < N), and is not an eigenvector 
of the transfer matrix in general. However, if the parameters Vj satisfy a set of constraints 
called the Bethe ansatz equation, the IV-particle state becomes an eigenstate. We call the 
state |^({u}iv)) on-shell state if {u} satisfies the Bethe ansatz equation, and off-shell state 
if no constraints are imposed on {«}. We then restrict the analysis to t = —1, which is 
equivalent to considering the six-vertex model (|2.14p , and prove the determinant form for the 
scalar products (^({u }iv) l*(M jv)) which is the inner product between the off-shell state 
|*(M*)) and the dual off-shell state (T({u}jv)|. 

From the L-operator, we construct the monodromy matrix 


M 

T a {v,{w}) = Y\£aj(,v/w M +l-j) 

3 =1 


(A{v,{w}) B(v,{w})\ 

\C(v,{w}) V(v,{w})J ’ 


(3.1) 


which acts on W a <S> (F\ <8> ■ ■ ■ <8> Fm)- Here we introduced the inhomogeneous parameters 
wi ,..., wm € C. See figure [2] for the graphical description of the elements of the monodromy 
matrix. 

Taking the homogeneous limit Wj —> 1 (1 < j < M), (12.121) is recovered: 


Ta(v, {ll’})|u,i = l,...,tt> M =l — Ta{v). 


(3.2) 


As in the above equation, hereafter we will omit {in} for the quantities in the homogeneous 
limit (e.g. A(v) := A(v, {h’})|k; 1 =i,...,iu m =i). Tracing out the auxiliary space, one defines the 
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transfer matrix r(e, {re}) € End(J r ® M ): 

t(v, {re}) = Tr Wa T a (v, {re}). (3.3) 

The repeated applications of the i?LL-relation leads to the intertwining relation 

Rabiu , v)Ta{u, { w})Tb(v , {re}) = %(v, {w})T a (u, { w})R ab {u , v). (3.4) 

Some of the elements of the intertwining relation are 
C(u, {w})B(v, {re}) - tB(v, {w})C(u, {re}) 

= g(u, v)(A(v, {w})V(u, {re}) - A(u, {w})V(v, {re})), (3.5) 

A(u, {w})B(v, {re}) = f(u, v)B(y, { w})A{u , {re}) + g(u, v)B(u, {re}).4(e, {re}), (3.6) 

V(u, {w})B(v, {re}) = f(v, u)B(v, {re})P(rr, {re}) - g(u, v)B(u, { w})V(u , {re}), (3.7) 

[B(u, {w}),B{v, {re})] = [ C(u , {re}), C(v, {re})] = 0, (3.8) 

where 

„. . ut — V , . (1 — t)v 

f( U ,v) = -, g(u,v) = -. (3.9) 

u — v u — V 

The transfer matrix r(e,{re}) is then expressed as elements of the monodromy matrix: 

r(e, {re}) = Triy a T a (v, {re}) = A(v, {re}) + V(v, {re}). (3.10) 

The arbitrary IV-particle state |^({e}jv)) (resp. its dual (\H({e}jv)|) (not normalized) 
with N spectral parameters {e}/v = {ei,... ,vjy} is constructed by a multiple action of B 
(resp. C) operator on the vacuum state |12) := |0 M ) := |0)i <S>• • -<g> |0 )m (resp. (f)| := (0 M | : = 

i(0|®---®Af(0|): 

N N 

I^(Mjv.M)) = JJ M)l n )> <^({^}rv, {rrr})| = (fi| JJ C(v j: {re}). (3.11) 

3 =1 3 = 1 

By the standard procedure of the algebraic Bethe ansatz, we have the followings. 

Proposition 3.1. The N-particle state |T({e}iv, {re})) and its dual (\H({e}jv, {re})| become 
an eigenstate (on-shell states) of the transfer matrix (13.1011 when the set of parameters {e}jv 
satisfies the Bethe ansatz equation: 


q(ej,{re}) 
d(vj, {re}) 
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f(vj,v k y 
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Then the eigenvalue of the transfer matrix is given by 

N N 

t(v, {re}) = a(v, {re}) JJ f(v, vj) + d(v, {re}) f(v jt v). 


(3.14) 
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For the case t = —1 which we investigate extensively in this paper, the Bethe ansatz 
equation (13.121) is 

M , \ 

II (— -P) =(- 1 ) iV+1 > j = (3.15) 

This is the Bethe ansatz equation for free fermions in the homogeneous limit Wj —> 1 (1 < j < 
N ). In the analysis below, we do not impose these constraints between the spectral parameters 
{u} and inhomogeneous parameters {re}. We remark that from this consideration on the 
Bethe ansatz equation, one can imagine the wavefunction is given by the Schur polynomials. 
The amazing thing is that a detailed analysis of its realization as partition functions lead us 
to a new combinatorial formula for the Schur polynomials itself. 

The scalar product between the arbitrary off-shell state vectors, which is mainly consid¬ 
ered in this section, is defined as 


(* (Mtv , M) I ^(Mtv , M)) 


N N 

(fl| C{uj, {w}) B(v k , {w})|ft) 

j=l k= 1 


(3.16) 


with Uj,Vk £ C. 

From now on, we specialize the parameter t of the f-deformed boson algebra to t = — 1. 
This is equivalent to considering the six-vertex model (12.141) . The following determinant 
formula in the homogeneous limit Wj —> 1 (1 < j < N) is valid. 

Theorem 3.2. The scalar product (13.161) in the homogeneous limit Wj —> 1 (1 < j < IVj is 
given by a determinant form: 


N 

<*c{«}jv)i*({»)jv)>= n<^) n 

j =1 1 <j<k<N 


(Uj +u k )(yj +v k ) 
(Uj - u k )(v k - Vj) 


detTvQUujivIMjv)) 


(3.17) 


where {u}n and {u}Ar are arbitrary sets of complex values (i.e. off-shell conditions), and Q 
is an N x N matrix with matrix elements 


V k - Uj 


(3.18) 


Here we will show the above determinant formula by using a method initiated by Izergin- 
Korepin EH m for the domain wall boundary partition function and recently developed 
by Wheeler [28] in the calculation of the scalar product of the spin-1/2 XXZ chain. This 
procedure was applied to a family of integrable five-vertex model [29] which, in contrast to 
the spin-1/2 XXZ chain, there was no need to impose the Bethe ansatz equation (i.e. on-shell 
condition) to show the determinant formula. For the six-vertex model we consider in this 
paper, we also do not have to impose the Bethe ansatz equation, i.e., the determinant formula 
(13.171) is valid for arbitrary off-shell states. 

First let us introduce the following intermediate scalar products which plays the key role 
for the proof 
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Figure 3: The graphical representation of the intermediate scalar products (13.191) with inho¬ 
mogeneous parameters WM+i-j associated with site j. 


See also figure 01 The term “intermediate” stems from the fact that (|3.19l) interpolates the 
scalar product (n = N ) (|3.16l) and the domain wall boundary partition function (n = 0) (13.281) 
(see also figure 01) . We have the following lemma regarding the properties of the intermediate 
scalar product. 

Lemma 3.3. The intermediate scalar product (13.191) <S , ({it}n|{^}jv|{ u; }) satisfies the follow¬ 
ing properties. 

1. S(Mn|M jv|{u>}) is symmetric with respect to the variables {u>i,..., wm-n+u} ■ 

& Y\f=M-N+n+l( l + P U n/™j)~ l S{{u} n \{v} Ar|{rc}) is a polynomial of degree M — N +n —1 
in u n . 

3. The following recursive relations between the intermediate scalar products hold 
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4■ The case n = 0 of the intermediate scalar products has the following form: 
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Proof. Property 1 follows from the LLL-relation 


Rjk(wM+l-j/ W M +l-k)L a k(u/ W M +l-k)Laj {u/WM+l-j) 

= L a j{u/wM+l-j)L a k{u/wM+l-k)Rjk(wM+l-j/w M+ i-k) 


holding in End (W a ® Vj (g> 14). Here R is given by 


R(v) 


(l 0 0 0\ 

0 f3(v — 1) v 0 

0 1 0 0’ 

\0 0 0 v) 


(3.22) 


(3.23) 


which intertwines the L-operators acting on a common auxiliary space (but acting on different 
quantum spaces). Note the usual RLL- relation (|3.22l) intertwines the L-operators acting on a 
same quantum space but acting on different auxiliary spaces. The above RLL- relation (I3.22|) 
allows one to construct the monodromy matrix as a product of the L-operators acting on the 
same quantum space (see also Appendix for an example of using its property to examine the 
symmetries of the domain wall boundary partition function), and rewriting the intermediate 
scalar products in terms of the resultant monodromy matrices makes one see Property 1 
holds. 

Property 2 can be shown by inserting the completeness relation into the intermediate 
scalar products 

n N 
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Property 3 can be obtained by setting u n = /3~ 1 wm-n+u in (I3.24|) . or can be directly 
observed by its graphical representation that the top row is completely frozen. 

To show Property 4, we first note by the graphical representation (see figure 0J) that 
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Figure 4: The intermediate scalar products (13.121 1) for n = 0, which corresponds to the domain 
wall boundary partition function. 


where 


N 

B N (v k ) = a ( 0 | ]^[ £aj(y k /w M +l-j)\l)a- (3.27) 

3 =1 

One can evaluate the domain wall boundary partition function (l iV | - =1 Bpf(vj) |0 jV ) by the 
standard procedure following the arguments of Izergin-Korepin M EZ], which is given in 
Appendix. The result has the following simple factorized form 
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(i K in^wi» N >=i 


2v 
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3 =1 


j =1 W M—N+j 1 <j<k<N 


3 


(Vj +V k ), 


(3.28) 


which together with (13.261) proves Property 4. □ 

Lemma 3.4. The properties in Lemma \3.3\ uniquely determine the intermediate scalar prod¬ 
uct (13.191) . 

Due to Lemma 13.41 the following determinant representation for the intermediate scalar 
product is valid. 

Theorem 3.5. The intermediate scalar product 5({u} n |{r;}7v|{u;}) (13.191) has the following 
determinant form: 


N _ 
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X TT Uj + Uk det N Q({u} n \{v} N \{w}) 
i<t<k<n U i~ Uk 


1 

Wk/Wj) 


(3.29) 
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with an N x N matrix Q({ti}n|{^}jv|{^}) whose matrix elements are given by 
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n R u - — w , 

l=M-N+n+l r 3 1 


/3uj + wi a(uj, {w})d(v k , {w}) - a(v k , {w})d(uj, {w}) 


V k - Uj 


1 


M 


WM-N+j 

l^M-N+j 


n (>-£'). 


, (l < j < n) 
(n+l<j<N) 
(3.30) 


Proof. We can directly see that the determinant formula (13.291) satisfies all the properties in 
Lemma 13.31 To show Property 4, we utilize the Cauchy determinant formula 


detjv 


1 


rii<j<fc<jv( a:: fc x j)(yj yk) 


- x 3~yk) 11 >/, 

Finally due to Lemma 13.41 the determinant formula (13.291) holds. 


(3.31) 

□ 


Corollary 3.6. Taking n = N in (I3.29P yields the determinant representation of the scalar 
product for the six-vertex model with inhomogeneous parameters (13.161 ) : 

(^({u} N , M)|^(Miv, {w})) 

=n^ n ^ d^q(w y iw,i W ) ( 3 . 3 2 ) 

3=1 1 <j<k<N v 3 K 33 

with 

0 (W»IW»IWU = ■faMUh.MI-lt.MHIi.M) (3.33) 

v k - Uj 

Further taking the homogeneous limit wj —> 1 (T < j < n) yields (13.171) in Theorem I S. .51 


4 Matrix elements of the t-deformed boson model and the 
six-vertex model 

In this section, we derive matrix elements for the t-deformed boson model for the generic 
parameter t, and then we restrict ourselves to the case of the six-vertex model t = — 1. 

Consider the arbitrary off-shell state, i.e., the parameters {u}tv in the IV-particle state 
(13.111) are arbitrary. The orthonormal basis of the IV-particle state |4 , ({u}at)) and its dual 
(T({u}at)| is given by \{n} M ,N ) ■= |ni)i u m )m and ({ u} m ,n\ -=i (ni[ <g>-• -® M («m|, 

where n\ + • • • + um = N. The wavefunctions can be expanded in this basis as 

l'I'(Wv)) = ^2 ({n}M,N\4’({v}N))\{n}M,N), (4.1) 

ni-|- \-n M =N 

(^(Miv)| = ^22 ({ n }M,N\{4>({v} N )\{n} M ,N}- (4.2) 

0<n\,...,n M <N 
ni -(- \-n M =N 
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There is a one-to-one correspondence between the set {h}m,n = {ni, • • •, um} {ni+- ■ ■+um = 
N) and the Young diagram x = (x\,X 2 , ■ ■ ■, xjv) (M > x\ > X 2 > • • • > xjsr > 1). Namely, 
each Young diagram x under the constraint x\ < M, i(x) = N can be labeled by a set of 
integers {n}M,N as x = (M nM ,..., l ni ). 

The following definition |23] on the ordering on the basis of particle configurations is 
useful for later purpose. 

Definition 4.1. [23] For two configurations {m}M, n+i = {mi, • • • >mjf} (mi + • • • + rriM = 

N + 1) and {n} M ,N = {«i, • • •, n M } {n i H-h n M = N), let Y1J = Y.k=j m k and YJj = 

J2k=j n k- We say that the particle configurations {m}M,N+i and {u}m,n are admissible, if 
and only if 0 < (ET - E") < 1 (1 < j < M ), and write this relation as {m}M,JV+i>{^}M,iV- 

Moreover we also define the ordering on the Young diagrams. 

Definition 4.2. For two Young diagrams y = (t/i, y 2 , ■ ■ ■, Vn+i) and x = (x\,X 2 , ■ ■ ■ ,xn), 
we say that y and x interlace, if and only if y 3 > Xj > Vj+i (j = 1,..., N), and write this 
relation as y y x. 

Proposition 4.3. Let {m}M,N+i and {u}m,n be the particle configurations described by the 
Young diagram y = (yi,..., yjv+i) and x = (x \,... ,X]\r). Then 

V y x {m} M ,N +1 > \n\ m,N• (4.3) 

For {m}M,N+i and {n}M,N, we introduce {p} r = {1 < p\ < ■ ■ ■ < p r < M} to be the 
set of all integers p such that m p = n p + 1, and {y} s = {1 < q\ < ■ ■ ■ < q s < M} to be 
the set of all integers q such that m q + 1 = n q . When {i7i}m,n+i and {u}m,n satisfy the 
admissible condition {m}M, n+i {p}r and {y} s satisfy s = r — 1 and Pk < Qk < Pk+i 

(k= l,...,r- 1). 

From the matrix elements of the L-operator, one finds 

{{m} M ,N+i\l3(v)\{n}M,N) = 0, unless {m} M ,N+i> (4.4) 

When the admissible condition is satisfied, one finds the following. 

M 

({ m}M t N+l\B{v)\{n}M,N) =a(0|({m}M,7V+l| £aj(v)\l) a \{n}M,N) 

3 = 1 

r r Qj — 1 

=t ,E'=i p J|(! _ TS-+ 1 ) n ll (1 -f3vt nk ), (4.5) 

3 =1 J = 1 k=pj +1 


where qo = 0, q r = M + 1. 


12 


This can be shown by combining the following partial actions: 


% 

-£aK' u )|l)as <S> | 

i = qj— 1+1 

= vPi-U-i(l-t n *i +1 ) JJ (1 - Pvt ni )\l) a {®l j =q ._ 1+1 \m k ) k } (l<j<r-l), 

i=Pj +1 
M 

C a l{v)\l) a ®{®kL qr _ 1+1 \n k ) k } 

l=q r — l+l 

M 

= _(%,+!) J| (1 - /3rt”0|0)„ {»,V ltl k)l} • (4.6) 

l=Pr +1 

Next, we examine the matrix elements of the one-row B and C operators furthermore 
at the point t = — 1. We first reduce the matrix elements (14.5 jl to a simpler form and then 
translate into the language of Young diagrams. The result for the matrix elements in the 
language of Young diagrams can be summarized as follows. 

Proposition 4.4. The matrix elements {{m}M,N+i\t3(v)\{n}M,N) at t = —1 are given by 


{{m}M,N+i\B{v)\{n} M ,N) 

1 + 2 p z f 2(1 + p z ) 

~ (1 + 0z) M + l \ 1 + 2 f3z 


#(y\ x ) 




N 

JX{1 + 2/3^(l - S Xj y j+ 1 )}. 


3 = 1 


(4.7) 


when x = (x\,... ,xn) = (M Um ,..., l ni ) C M N and y = (yi ,..., vn+i) = (M mM , ..., l mi ) C 
are strict partitions satisfying y + x, and zero otherwise. Here jf(y\x) denotes the 
number of parts in y which are not in x. We regard the product in the right hand side of 
(021) as 1 when x = <fi. 

Proposition 4.5. The matrix elements ({w} M,N\C(v)\{m} m,n+i) at t = — 1 are given by 


({n}M,N\ C ( v )\{ m }M,N+l} 

1 \ 2(1 + P z ) 

z{l + /3z) M - L \ 1 + 2 fiz 


#(w v l* v )- i 


y^iV+1 

z ^j=l 


, , v _W JV r ' 
Vi 2^,=i x - 


N 

H{l + 2Pz(l-6 x y y v +i )}. (4.8) 

3 =1 


when x v = (M + l^/x = (x^,... , x)y) = (M ni ,..., l njVf ) c and y v = (. M + l) N+1 /y = 
(yf ,..., y)v+i) = (-^d m S • ■ •, l mM ) C M w+1 are sfricf partitions satisfying y v + x v , and zero 
otherwise. We regard the product in the right hand side of (14.81) as 1 when x v = <f. 

Proof. We show 021). Eq. (14.81) can be calculated in the same way. 

First, we note that since n/= i(1 — t nv i +1 ) / 0 unless n p . = 1 (mod 2) for all j, we 
only need to consider the case when both {m}M,N +1 and {u}m,n are sequences of numbers 
0 and 1. Otherwise, the matrix elements vanish as explained in section 2. Translating this 
restriction to the language of Young diagrams, this means that we restrict both the partitions 
y and x to be strict. 
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The matrix elements for the case y >- x is calculated as (14.51) for generic t, which can be 
furthermore simplified at t = — 1 as follows 

({m}M,N + i\B(v)\{n} M ,N) =2 r «S=i«-^«(l - / 3 u )#{ fceu J=ite+ 1 --,®-i}l^=o} 

x (1 + 1 }l”fc= 1 l ) ( 4 . 9 ) 

where {m}M, jv+i and {u}m,n are both sequences of 0 and 1 satisfying {i7i}m,n +i > {n}M,iV- 
Using 

(1 — p v ^#{ kGU j=i{Pj +1 --’ c U- 1 }\ n k=°} 

— (1 —/3u) # ^ eL J = lte + 1 >-.®- 1 l}(l - / g t ,)-#{ fceU J =1 {p j + l,...,qj-l}|n fc = l} 

= (1 — /3 W )££i ^—SJ=iPj+ M + 1 -r-#{feeUj =1 {p ;l +l,...,g i -l}|n fc =l}^ (4.10) 

({m} M ,iV+i |£(v) | {n} M ,iv) 

=2 r v EJ=i Pi-E* (1 _ ^)E^i ®-EJ=i W+M+l-r 


one has 


/ ^ _|_ Qy \ #{fc€ur_i{pj-+i,...,gj i}|nfc—1} 

x VT^J 

=2 r (l - i 5u) M+1 - r (u- 1 - /3)££ifc-Ei=i« 

' ^ -|- j$y \ l}| n fc = 1} 


-/w 


(4.11) 


From the translation rule m 

r r— 1 7V+1 iV 

E^-E v = E - E 

1=1 1=1 /—i 1=1 

r - 1 + #{& € U^ =1 {pj + 1, ■ ■ ■, qj - l}\n k / 0} = #{j G {1,..., -/V}^ / %+i}, (4.12) 

One gets 

({"i}m,at+i |#(v) |{h}m,at) 


=2 r (l - ( Su) M - r+1 (u- 1 - /3)Ef= 


, ,,- E ^ % / 1 + M #Oe{b-^}|^M + i}-r + l 

Vi-^/ 


=2 r 


l + /3u 
1 -/3u 


1 —r 


(1 - Pv) M+1 ~ r {v~ l - / 3 )E J= 1 ^-E j= i % 


TV 




1=1 


(4.13) 


Introducing the variable z = (y 1 — /?) 1 and changing the variable from v to z, the matrix 
elements can be rewritten as 


{{m}M,N+i\B{y)\{n}M,N) 

= 2 r (l + 2Pz) 1 ~ r (l + Pz) r - M ~ L z^ 


1 —r (i , 


AT 


j = l ^ 


nii+wi-wi. ^- i4 ) 

1=1 
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where r =: ff{i)\x) is the number of parts in y which are not in x. 


□ 


Note that the elements of the T-operator f-'ajiy) m (ESI) at t = — 1 reduces to those 
for the six-vertex model (|2.141 ) : 


C a j{v) = 


(l-Pv 

0 

0 

V o 


0 


0 


1 + f3v 2v 


v 

0 


°\ 

0 

0 

v) 


aJ 


1 +/3z 

0 

0 

V 0 


0 

1+2/32 

l+0z 

1 

0 


0 

2 2 

1 +0Z 

z 

1 +0Z 

0 


0 \ 
0 
0 

T Tpl) 


(4.15) 




Example 4.6. We set M = 5 and N = 2 and consider the case {m}M,N+i = {10011} and 
{n}M,N = {00110}. Translating into the language of Young diagrams, we have y = (5,4,1) 
and x = (4,3). One finds ff{y\x) = 2 since 5 and 1 are in y but not in x. We also have 
Y^=iVj ~ J2j=i x j = 3 and x\ = j/ 2 > x 2 / V 3 - The right hand side of (14.7|) becomes 
2 2 z 3 (1 + /3z)~ 4 which can be easily checked to match with the left hand side. 

For /3 = 0, the matrix elements reduce to the skew Schur Q-polynomials [241 |2S] 

({m}M,N + i\B(v)\{n} M:N ) = 2#^z^=" y ^U*i = Q y/x {z), (4.16) 

which we will use to derive the wavefunction in the next section. 


5 Wavefunctions of the six-vertex model 

Let us examine the wavefunction at the point t = — 1 where the f-boson model reduces 
to the six-vertex model. We show the corresponding wavefunction is essentially the Schur 
polynomials. 

Definition 5.1. The Schur polynomial is defined to be the following determinant: 

det N (z* k+N ~ k ) 


s\(z) = 


(5.1) 


ril <j<k<N( z j z k) 

where z = {zi ,..., zn} is a set of variables and A denotes a Young diagram A = (Ai, A 2 ,..., Aat) 
with weakly decreasing non-negative integers Ai > A 2 > • • • > Xn > 0 . 

We show the following equivalence between the wavefunction of the six-vertex model and 
the Schur polynomials. 

Theorem 5.2. The wavefunction has the following form 


N 


2 W nJ=1 Z j Ell <j<k< N ( Z j + Z k + 2 fiZjZk) 


{{m} M ,N | [ 
3 = 1 


nf=i(i + 13*3 




s\(z). (5.2) 


where Zj = (■ v ■ — /3) 1 and A = (Ai, ... ,Xn) C (M — N) n is a partition related to x = 


.— 1 ,3 )~ 1 

(xi, .. .,x N ) = (M mM ,... ,l mi ) by A j = Xj - N + j - 1. 
The dual wavefunction has the following form 
N 


/r ,i TT nt nr l \ Hl<j<k<N( z 3 + z k + 2f3zjZ k ) 

wIUC^OIWm.jv) = — = n Ai - i _ v -s A v(«), 


j =1 Tlj=l0- + P z j) 

where A v is the Poincare dual of A: A v = (A{,..., A)y), A J = M — N — A/v+i-j. 


(5.3) 
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We remark that the integrable model we consider here seems to be a special case of the 
one considered in [12] , whose wavefunction was obtained by the coordinate Bethe ansatz up to 
a normalization factor. However, to prove combinatorial formulae, it is crucial to determine 
the exact form starting from the first principle, i.e., starting from the L-operator, since we 
combine the above theorem with the exact expression for the matrix elements (14.71) and (14.81) 
derived in the previous section to derive a new combinatorial formula for example. 


Proof. We show (15.21) . Eq. (15.31) can be proved in the same way. First, we redefine the 
L-operator as 


C aj (z) = (1 + f3z)£ aj (v) 


(l 0 0 0 \ 

0 1 + 2/3 z 2z 0 

0 1 + Pz z 0 

\0 0 0 z) 


and the corresponding monodromy matrix 

M ~ (j 

7a( Z ) = n^W = 

and show the following equivalent equality for (15.21) 

N N 


B(z)\ 

V(z)J 


{{m}M,N | n+,)in>=2«n^ n (Zj + z k + 2f3zjZk)s\(z). 

j =1 j =1 1 <j<k<N 

To prove this, we first show the following lemma 


(5.4) 


(5.5) 


(5.6) 


Lemma 5.3. 

Ui< j<k < N ( z : + z k + 2/3 ZjZk) 

does not depend on /3. 


Proof. We prove this lemma by showing the following properties for ({m}M,N | n)=i B{zj)\£l)\ 

1. ({m}M,N | nf=i &( z j) 1^) is a polynomial of /3 with highest degree N(N — l)/2. 

2. ({m}M,N | nf=i &( z j) 1^) has Zj + Zk + 2 /3zjZ k , 1 < j < k < N as factors. 

We first show degp({m}M,N\ nf=i £>(zj)|fi) — 3V(iV — l)/2 by induction on N. The case 
N = 1 follows as an special case of the general fact degp({m} M,N+i\B(z)\{n} m,n) < 3V 
which can be seen easily from the definition of the L-operator £(z). Next, let us assume 
deg p({m}M,N\Uf =1 B( Zj m < N(N-l)/2. One can see deg^ {{m} M , jv+i| II f=i &( z j) 1^) < 
(N + l)3V/2 by combining the assumption deg^ ({m}M, n\ n,=i &i z j)\ty — N(N — l)/2, the 
fact deg^{{m} M,N+i\B{z)\{n} m,n) < 3V and the decomposition 


N +1 

({m} M ,N+ 1 | Q B{zj)\Vt) 


3=1 


N 

^ ({m} M ,N+i\B(z N+1 )\{n} M ,N){{n}M,N\ £>(+,)| 0 ). 

{"■1m, n i=i 


(5.8) 
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Next, we show Property 2. It is enough to show the case N = 2. The case for generic N 
follows from the commutativity [B(zj),B(zk)] = 0 and the decomposition 

N N 

{{m} M ,N | ^ ({m} M ,N\ni3( z j)\{n}M,2)({n} M } 2\B( z i)B( z 2)\ty. (5.9) 

j = 1 Wm, 2 f=3 

Now we show for the case N = 2 by induction on M. Let us denote the A, B, C, D op¬ 
erators consisting of M L-operators as E>m(z) for example. The case M = 2 can be 
checked explicitly {1,1\Bm{zi)Bm{,z 2 )\Q) = 4zi z 2 (z\ + z 2 + 2(3z\z 2 ) ■ Let us assume that 
(xi, X 2 \Bm(zi)Bm(z 2 )\&) has zi+Z 2 + 2 / 3 ^i 22 asafactor. We examine {x\, X2 \Bm+i(zi)Bm+i{z 2 )\£1)■ 
One can easily show by its graphical description that 

{xi,X2\Bm+i{zi)Bm+i{z 2)\®) = {ziZ2) Xl ~ 1 (xi,X2\B X2 - xl+1 (zi)B X2 - xl+1 ( y Z2)\V). (5.10) 

If xi / 1 or X 2 / M + 1, {x\,X 2 \Bm+i{zi)Bm+i(z 2 )\Q) has z\ + Z 2 + 2f}z\Z2 as a factor by 
assumption. 

We examine the remaining case x\ = 1,xm+i = M + 1. We show 

~M _ „M 

(1 ,M + l\B M +i(zi)B M +i(z 2 )\ty = / Lz 1 z 2 — - —{zi + z 2 + 2/3ziz 2 ). (5.11) 

Z\ — Z 2 

By graphical description, one sees 

{xi,X 2 \B M +i(zi)B M +i(z 2 )\£l) = 2zi{f M (zi,z 2 ) + 2z| f+1 (l + 2/3zi)}, (5.12) 

where /m( 2 i, 22 ) = (1|Dm ( 2 i)jB]vf ( 22 )!^). Again, we use its graphical representation to derive 
the following recursive relation 

fM(zi,z 2 ) = zi{f M -i(zi,z 2 ) + 4^(1 + /Lzi)}, (5.13) 

with the initial condition 

/2 ( 21 , 22 ) = 2^iz 2 (2i + 2z 2 + 2/32Ti^ 2 )- (5.14) 

We can show by induction that 

M _ M 

/m( 2 i, 2 2 ) = 2z 2 —- —(21 + 2 2 + 2/Lziz 2 ) - 2z^ +1 (l + 2/3zi), (5.15) 

2 l — 22 

solves the recursive relation (|5.13|) and the initial condition (15.141) . Hence, the expression 

(15.111) follows from (15.121) and (15.151) . We thus have shown by induction that 

(x±, x 2 \Bm+i(zi)Bm+i(z 2 )\£1) has 21 + z 2 + 2f3ziz 2 as a factor, and Property 2 is proved. 

From Property 2 we have, deg^({m} m.n\ Il/Li B( z j)\Q) > N(N — l)/2. Together with 
deg^({?7i}Af,iv| UU < 1V(AT —1)/2 which we proved before, we have Property 1. □ 

From Lemma 15.31 one can examine 

(W Af ,iyinjL 1 g( 2 j)|n) l 

W.i<j<k<N^ z i 4" z k + 2 fdz.jZk) 

by dealing the case j3 = 0. 
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(5.17) 


Lemma 5.4. We have 

({m} M ,N injii*?(zj)i n ) 

rii<i<fe<Ar( 2 j "t Zk "t 2(3 z jZk) 

Proof. To prove the lemma is equivalent to show 
N N 

{{m} M ,N\ H#(zj)|fi)b=o = 2N Y{ z i n ( z j + Zk)s\{z) = Q x (z). 

j =1 j =1 1 <j<k<N 

where Q x ( z ) is the Schur Q-function 

Qx(z) = 2 n ^2 w 

u!S5jv 


N 


riv n 


, z j 

ij =1 l<j<fe<iV 


Zj + Zfc 
z? - Zk 


N 

^nww). 


(5.18) 


(5.19) 


where 2 = {zi,..., 2 jv} is a set of variables and x denotes a strict Young diagram x = 
(xi, X 2 , ■ ■ ■, xn) with strictly decreasing non-negative integers x\ > X 2 > ■ ■ ■ > xn > 0 . 
( 15 . 181 ) follows from the fact ( 14 . 161 ) derived in the previous section that the matrix element of 
a single £>-operator at /3 = 0 is nothing but the skew Schur Q-function 


{{m}M,N+l\B(z)\{n}M,N)\p=0 = ({"l}M,7V+l|l?('y)|{^}M,7v)|/3=0 

= 2 #(vl*) z EiL + 1 1 Wi-EjLi*j 


= Q y/x (z). (5.20) 

(15.181) follows as a consequence of the addition formula for the skew Schur Q-function 

{{m} M ,N+2\l3(z 1 )B(z2)\{n}M,N)\i3=o 

= y^ J {{m} M ,N+2\B(zi)\{£} M ,N+l)\/3=o{{P}M,N+l\B{z2)\{n}M,N)\l3=0 

e 

^ ^ Qy/w(z\ )Q W / X {Z2) 

W 

Qy/x{zi, z 2 )• (5.21) 

□ 

Combining Lemma l54l and (15.171) . we have (15.61) . and the proof of (15.21) is completed. □ 


6 Combinatorial formulae for the Schur polynomials 


By combining the analysis of the partition functions in the previous sections, we obtain 
combinatorial formulae for the Schur polynomials. 

Theorem 6.1. We have the following combinatorial formula for the Schur polynomials 

s ^ z ) =77— 1 .., o Q - „ ^ e n 




111 <j<k<N( z i + Zk + 2/3zjZ k ) 


2(1 + (3zk)\ 

1 + 2/3 Zk ) 


x( N )yx( N ~ •*■)>— yx(®)=(j> k—1 
#(a;( fc )|a:( fc - 1 ))-l k -1 


JJ[ ( 1 + 2/3zfc(l — h^k- 

3 = 1 


i) r W , 


( 6 . 1 ) 
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where (3 is an arbitrary parameter, x = (x^\ ..., x^), k = 0,1,..., N are strict partitions 
satisfying the interlacing relations x^ N ^ + x^ N ^^ >-■■■>- x^ = 4>, and x^ is fixed by the 
Young diagram A = (Ai,... , Ajv) as \j = Xj N ^ — N + j — 1 . 

Proof. We decompose the wavefunction as 


({™}M,7v|'i'(WA0) 

N 


We insert the expression (15.21) in the left hand side of (16.21) on one hand. On the other hand, 
the right hand side of (16.21) can be expressed using the evaluation of matrix elements (14.71) . 
Combining the two expressions and simplifying gives the combinatorial expression for the 
Schur polynomials (16.11) . □ 


Example 6.2. Let us check the case N = 2, A = (1,0). x^ is fixed as x (2) = ( 1 , 0 )+ ( 2 , 1 ) = 
(3,1) and x ^ = </>. a^ 1 ) satisfying the interlacing relation x^ + a^ 1 ) + a:^ 0 ^ has three cases 
x W = (1), a;*- 1 ) = (2) and a^ 1 * = (3). Each term in the sum of the right hand side of (16.11) 
has the contribution 


z 


z 


z 


4-1-1 

2 


4-2-1 

2 


4-3-1 

2 


2(1 + fiZ2) N 

l+2fiz 2 y 

1-1 

) Z 1~ X , 

(6.3) 

2(1 + fizfi ) N 

1 + 28z 2 y 

2-1 

) (1 + 2 / 3 ^ 2)21 1 ) 

(6.4) 

2(1 + fizfi ) N 

l+2fiz 2 y 

| (1 + 2fizfi) z\ 1 , 

(6.5) 


which sums up to (z\ + 22 + 2fiz\zfi){z\ + zfi). Dividing by z\ + 22 + 2/3ziZ2, we have z\ + Z 2 , 
which is nothing but S(i,o)(- 2 i, 22 ). 

Proposition 6.3. The following well-known identity holds true for the Schur polynomials. 

1 


Y s x (z)s x v(y) = 


\C(M-N) N 


l<j<k<N 


(Zk - Zj){yj - y k ) 


det^r 


„M __ M 
~j 

Zj - Vk 


( 6 . 6 ) 


Proof. First, substituting the completeness relation, one decomposes the scalar product as 


(lk({u}Ar)|4'({D}Ar)) = ^ WMaOI Wm,7v)({™}m, 7V ^(Wiv))- (6-7) 

Then substituting the determinant representation for the scalar product (13.171) into the RHS 
of the above and utilizing the relations (15.21) and (15.31) yields the Cauchy identity (16.61) by 
changing the variables from Uj and Vj to yj = ( uj 1 —/3) _1 and Zj = ( vj 1 —/3)" 1 respectively. 

□ 
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7 Enumeration of alternating sign matrices 

In this section, we make an application of the domain wall boundary partition function to 
the enumeration of alternating sign matrices. See [3ol EU [321 E3 EH ESI ESI E3 [381 [39] for 
example of a huge literature on the relation between the enumeration of alternating sign 
matrices and integrable vertex models. The presentation below for the explanation of the 
relation between alternating sign matrices and the six-vertex model follows the lines of 135] . 
We take the homogeneous limit of the domain wall boundary partition function (13.281) or 

EB 


n n 

^(Mni{iD=(i n in^(^)i Q )=n^) n (^+^)- ( 7 -i) 

k =1 j =1 l<j<k<n 

An algebraic analytic proof of (13.281) is given in Appendix. It can also be obtained as a 
special case of the wavefunction (15.21) . 

As a corollary of the domain wall boundary partition function, we derive a simple expres¬ 
sion for a special case of the generating function of alternating sign matrices. 

Definition 7.1. Alternating sign matrices are square matrices with the following properties: 

1. each entry is either 0, 1 or — 1. 

2. there is at least one nonzero entry in each row and column. 

3. the entries in each row and column sum to 1. 


To satisfy the above conditions, the nonzero entries must alternate in sign along each row 
and column. 


Example 7.2. For n = 3, there are 7 alternating sign matrices: 

( (\ 0 0 \ /0 0 i\ (\ 0 0 \ 

ASM(3) = < 0 1 0,0 1 0 , 0 0 1, 

l\0 0 1/ V 1 0 0/ \0 1 0/ 

/0 0 1 \ /0 1 0 \ /O 1 0 \ /0 1 0 \ \ 

1 0 0,1 0 0,0 0 1,1 -1 1 

\0 10/ \0 0 l) \l 0 0/ \0 1 0/ J 


(7.2) 


Definition 7.3. Let ASM(n) be the set consisting of all n x n alternating sign matrices. For 
A G ASM(n), let us define 


v(A) := AijAht, (7.3) 

l<i<fc<n 

1 <£<j<n 


fj.(A) as the number of —l’s in A and p(A) as the number of 0’s to left of 1 in the first row 
of A. The generating function of the alternating sign matrices is defined as 


Z ASM (n,x,y,z) = ^2 x" (A V (A) £ p(A) - (7.4) 

AeASM(n) 


20 








Each element of an alternating sign matrix ASM(n) has one-to-one correspondence with 
a particular configuration which makes a non-zero contribution to the domain wall boundary 
partition function of the six-vertex model. To simplify the explanation, we identify a par¬ 
ticular configuration of the domain wall boundary partition function at the coordinate [j. k) 
with the matrix elements a(*|j(*|Ay(vfc)|*) 0 |*)j. 

Definition 7.4. Let 6VDW(n) be the set of all particular configurations making non-zero 
contributions to the domain wall boundary partition function of the six-vertex model. For 
C € 6VDW(n), Let i /(C) be the number of vertex configurations a (l| J (l|>C a j('yfc)|l) a |l)j, 
j,k = 1,... ,n in C, p(C) the number of vertex configurations a (l| ; (0|£ aj (ufc)|0) a |l)j, j , k = 
1 ,..., n in C, p(C) the number of vertex configurations a {l\j{l\C a j{y n )\V) a \^)j, j = 1,... ,n 
in C. 

Theorem 7.5. fJOty There is a bijection between the set of alternating sign matrices {A € 
ASM(ro) | v(A) = p,p(A) = m, p(A) = k} and the set of configurations making non-zero 
contributions to the domain wall boundary partition function of the six-vertex model {C € 
6VDW(n) | zz(C) = p,p(C) = m, p(C) = &} by identifying the matrix elements 1 with 
a (0|j(l|£ a:/ '(wfc)|l) a |0)j, —1 with a (l\j(0\£ a j(vk)\0) a \l)j and 0 with one of the rest of the four 
non-zero configurations. 

Proposition 7.6. fJUj Let #a{*\j{*\£aj(vk)\*)a\*)j be the number of vertex configurations 
a{*\j{*\£aj(vk)\*)a\*)j, j,k = 1,... ,n in a configuration C £ 6VDW(n). 


# a <0| J -(0|£ oj (w fc )|0) a |0) j =u{C), 

(7.5) 

# a (0| i (l|T a j(ufc)|0) a |l) j = n(n - l)/2 - u(C) - p{C), 

(7.6) 

#a(0| i (l|£aj(u fc )|l)a|0)j = p{C) + U, 

(7.7) 

# o (l| J -(0|£ < y(« fc )|0)a|l)j =P(C), 

(7.8) 

#a(l|j(0|£ O j(^fc)|l)o|0)j = n(n - l)/2 - v{C) - //(C), 

(7.9) 

# 0 (l| J -(l|£ <y -(« fc )|l)a|l) J - = v{C). 

(7.10) 


We use CCD, Theorem 1731 and Proposition 17.61 to show the following simple formula for 
a particular type of the generating function of alternating sign matrices: 

Proposition 7.7. 

Z ASM (n,u-hu,l) = u n ^-^ 2 . (7.11) 

Proof. We first take the homogeneous limit Vj —> v, j = 1,..., n of the partition function of 
the six-vertex model (17.111 : 

Z(v) Z(M„|{1})| m ,...,^ = (2i,)"<"+W 2 . (7.12) 

On the other hand, by the definition of domain wall boundary partition function and Propo¬ 
sition 1731 we have 

Z(v) = (aoaiY^^boh^-^^-^-^c^c^ 0 ^, (7.13) 

C€6VDW(n) 
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where ao = 1 — /3v, a\ = v, bo = 1 + /3v, b\ = v, cq = 1, c\ = 2v. 
Combining (17.121) and (17.131) and simplifying, one gets 



(7.14) 


We make change variable from v to u := 2/(1 +/3v). Finally, we use Theorem l7.5l on the one-to- 
one correspondence between the set of alternating sign matrices and the set of configurations 
of the domain wall boundary partition function to replace the sum over C € 6VDW(n) and 
the numbers i/(C), /i( C ) by the sum over A € ASM(n) and the numbers v(A), Then 

(17.141) can be rewritten as (17.111) . which concludes the proof. □ 

Example 7.8. For n = 3, there are 7 alternating sign matrices as in (17.21) . Summing up all the 
corresponding factors, we have Zasm( 3, u—1, u, 1) = (it—l) 3 +2(u— l) 2 +(it—l)it+2(u—1)+1 = 


We finally remark that it may also be possible to derive (17.111) as a limit of the determinant 
representation for Zasm(M, x, y, 1) in [38] . 

8 Conclusion 

In this paper, we derived a new combinatorial formula for the Schur polynomials. The quan¬ 
tum integrability is useful to derive a formula even for Schur polynomials which is the most 
fundamental symmetric polynomials. The formula expresses Schur polynomials with an addi¬ 
tional parameter besides the spectral parameter. The other known formula in a similar sense 
is the Tokuyama formula, which can be interpreted as a deformation of the Weyl character 
formula and the determinant expression. The representation theoretic meaning of the defor¬ 
mation parameter in our formula is unknown now which may be worth investigating. There 
may be other possibilities to find combinatorial formulae to express Schur polynomials and 
other symmetric polynomials in terms of additional parameters using the power of quantum 
integrability. This may be achieved by dealing with partition functions consisting of different 
local L-operators or changing global boundary conditions for example. 

Besides the traditional problem of the application of partition functions of the six-vertex 
model to the enumeration of alternating sign matrices, one of the most active line of researches 
on quantum integrable combinatorics today is to derive combinatorial formulae for symmetric 
polynomials such as the Cauchy identity, Littlewood identity and so on by analyzing the 
transfer matrices or partition functions of integrable lattice models. The power of quantum 
integrable combinatorics is that one can fuse combinatorics with algebraic analysis to find 
and prove identities which seems to be hard to show in a purely combinatorial or a purely 
algebraic way. See 0 El 1221 mi M for finite lattice and [33] for infinite lattice for example 
of the recent progresses on this line. 
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A Evaluation of the domain wall boundary partition function 


Here we derive the factorized form (13.281) of the domain wall boundary partition function 
defined by 

N 

Z({v} N \{w} N ) := <1*| [] B N ( Vj )\0 N ), (A.l) 

3 = 1 

where is given by (13.271) and {xr ;}tv = {u>m+i-n, • • • ,%}. 

Proposition A.l. The domain wall boundary partition function dS3D is expressed as the 
following factorized form: 

N 2v 

z{{ v }n\{w } N )=n —f n ( v j +v k)- A- 2 ) 

j =1 w j l<j<k<N 

Proof. This can be shown in the standard approach due to Izergin and Korepin [Ml EH- 
First, one shows the following four conditions. 

Lemma A. 2. The domain wall boundary partition function dSUD satisfies the following 
properties. 

1. YY',Lm+i-n w ^ + ^ AI Z({v}n\{w}n) is symmetric with respect to the variables {ic}at. 

2. Z({v}n\{w}n)/vn is a polynomial of degree N — 1 in vn- 

3. The following recursive relations between the domain wall boundary partition functions 
hold: 


Z (Mat | {w}n) | vn^-'wm-n+i 

M 

n 

j=M-N +2 


2 

1 


WM-N +1 
/ 3wj 


N -1 

n 

7 — 1 


0V] ) ^({uU-ilMiV-i). (A.3) 

wm-n+i J 


4- The case N = 1 of the domain wall boundary partition function following form: 


z(Mi|Mi) 


2v\ 

wm' 


(A.4) 


Properties 2, 3 and 4 can be shown easily with the help of the graphical representation (see 
figure El for Property 3) of the domain wall boundary partition function Let us explain Prop¬ 
erty 1. First note that the domain wall boundary partition function Z({v }n IM n) can be re¬ 
expressed using the transfer matrix Tj(wM+i-j) = fl^Li B a j(v a /wM+i-j) £ End(VF® JV <g) Tj) 
propagating in the horizontal direction as 


Z ({ v }n\{w}n) = a(l A \C(wm+1~n) • • • C(wm) |0 A )a, 


where 


Tj(wM+l-j) 


{A(w M +i~j) B(w M +i-j)\ 
\C(w M +i-j) D(wM+i-j)J ■ ’ 


(A.5) 


(A.6) 
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Figure 5: A graphical description of (1A.3D . 


From the RLL relation, one has 

C(w k )wjC(wj ) = C(wj)w k C(w k ). (A.7) 

Combining (1A.5D and (IA.7D shows Property 1. 

The remaining thing to do is to find the explicit forms of the functions satisfying the 
properties in the Lemma. One can easily show that 

N M 

Z{{v} N \{w} N ) = n<>i) II wf~ N ~ 3 JJ (Vj + V k ), (A.8) 

j =1 j=M+l—N l<j<k<N 

satisfies the above four properties. □ 
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